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Abstract. Some questions concerning the calculation of the number of "physical" (metastable) states or 
complexity of the spherical p-spin spin glass model are reviewed and examined further. Particular attention 
is focused on the general calculation procedure which is discussed step-by-step. 

PACS. 75.10.Nr Spin-glass and other random models - 02.30. Mv bho 

1 introduction of physical problems. In this context classical calculations 

pin?) for the complexity of the Sherrington-Kirkpatrick 0] 
The analysis of the equilibrium and non-equilibrium prop- (g^) and other disordered spin models have been recon- 



o 

O . erties m terms of the energy landscape origmally pushed gidered, extended and in some cases criticised 
. ^ , forward for the structural glass transition H , has risen in 

X 

^ [ the recent years a new mterest on the topological proper- Motivated by these criticisms in this paper some ques- 

ties of the energy or free-energy landscapes of disordered tions concerning the calculation of the complexity of disor- 

and complex systems. In this approach an important the- dered spin systems are reviewed and examined further in 

oretical tool is the logarithm of the number metastable a solvable model, the spherical p-spin spin glass (pSP-SG) 

states, called complexity or configurational entropy, iden- model introduced by Crisanti and Sommers 0. The aim 

tified as basins or valleys on the landscape. of the paper is not the calculation of the complexity for 

Solvable models, such as mean-field models, have al- the spherical pSP-SG model, which has been computed in 

ways played an important role in the theoretical study Refs. |??lll()ll29j . but the procedure of calculation itself to 



2 A. Crisanti et al.: The Complexity of the Spherical p-spin spin glass model, revisited 

clarify know results which could be useful for the under- the minimum (and its basin of attraction) is a metastable 

standing of the complexity of other systems. The spherical state of infinite life-time. However, despite this simple in- 

pSP-SG model is only used to enlighten subtle points of tuitive picture, not all minima of -F'tap('tT') can be asso- 

the procedure. We have tried to make the paper almost ciated with physical metastable states but only those for 

self-contained so that it can be also used by readers in- which (Plefka's criterion) |12II13| 

terested into the problem but not too familiar with all r(n) „ ^ 

:rp = l-^^(l-m,2)2>0 (2) 

reported results. 

™ 1 1 ; ,, , , ;• r 1 •; wher e /3 = 1 /T aud c ( (7 ) is a f uuc t iou of (7 = ( 1 / ) ■ TO^ 

ine general approach to the calculation of complexity -i \ , , 

. , i-n-ra whose form depends on the interactions. For example for 

m mean-held spin glass model is discussed m Section |2J 

„i , • -1, , ; , • ; • m • ;i 1 thc SK modcl c(g) = [T^ltT^ . while for n-spin inter- 

ihe procedure is illustrated m Section |4| using the spher- ""^ — ' 

• 1 QT3 cz-i J 1 • 4- J J • c r action models c{q) = (/3^p/2)(p - l)q^^^ HH. Therefore 

ical pSP-SG model introduced m Section |3| The of cor- ^ ' ^ ' ' ' ' 

, ... 1 • , • in n within this approach the calculation of the number of 

rectness the procedure is discussed m Section |5| l^mally 

„ . ™ , • 1 • 11- • metastable states is reduced to that of counting the num- 

Section |6| contains some conclusions and discussion. 

ber of solutions of J^l which are minima and satisfies the 
Plefka's criterion Q ( "physical" minima) . If it were found 

2 How to compute the Complexity 

that physical minima must satisfy additional constraints, 
The metastable states in mean-fields spin models are more those must be also included. 

easily studied using the Thouless- Anderson-Palmer (TAP) Different physical minima may have different free-energy 

method .11 , which introduces a mean- field free energy density, thus to have a better description of metastable 

functional Ftap("T') of the local magnetizations m = (rni, states one can group together all minima with the same 

7712, . . . ,TOAr), where N is the number of spins. At any free-energy density and introduce the function p(/) which 

temperature T the metastable states of the system are gives the number of metastable states with ^tap("t.) = 

identified with the local minima of Ftap("i), i.e., with Nf. The configurational entropy is then defined as: 

the solutions of i 

E{f) = -\np[f). (3) 

9™,i^TAp(m) = 0, i = l,...,7V (1) 

We are eventually interested into the large N limit, thus 
with the additional requirement that aZZ eigenvalues of the S{f) is different from zero only if the number of physical 
matrix 9™ . (9^ . -Ftap {'m) evaluated on the solution are pos- minima with free energy density / is exponentially large 
itive. At the mean-field level different local minima are with N . 

separated by infinite barriers therefore the system cannot If we label the Nso\ solutions of the TAP equation ^ 
escape from a local minimum in a finite time and hence with the subscript a {a = \, . . . ,7Vsoi) by definition p{f) 



A. Crisanti et al.: The Complexity of the Spherical p-spin spin glass model, revisited 3 
is given by same set of solutions, an assumption that must be verified 



a=l i=l 



in each case (and for each value of /) separately. 

Keeping the sum over all solutions is not, however, 
completely free of difficulties: since all solutions are counted 
where Af (i = 1, . . . , iV) are the eigenvalues of the Hessian the determinant of the Jacobian can be negative and the 



1 ; 1 L'' 

x9{x^)S[FTAp{m)-Nf] (4) 



matrix for the a-th solution: absolute value must be retained making the subsequent 

calculation more problematic. To overcome this difhculties 



Xfj = dm^dm, -PrAp(m)|^^^„ (5) 



the absolute value is simply dropped leading to expression: 



and 9{x) is the Heaviside theta-function. As it stands l@J ^ 

is difhcult to handle, however using the properties of deha- M) = 11 / dm, S [9™i-FTAP (m)] 

i—l 

function it can be transformed into the more manageable x det(x(m)) S [FTAp('m.) — Nf] , (7) 



form: 



N r 

p{f) = n [y dm,0(AO5[9„,FTAp(m)] 



and arguments are given to justify under which circum- 
stances this reproduces the correct result with the absolute 



1=1 

xdet(x(m))^^(xp),5[i^TAp(m)-iV/] (6) ^^1^^. 

where , the Hessian matrix © evaluated for a generic We have to compare ~p{f) given by with p{f). The 

m, is the Jacobian of the transformation and A, are its question is when ~p{f) yields the same resuh as p{f). The 

eigenvalues. The theta-functions ensure that the determi- ™ain difference between © and Q is the support of the 

nant of x is always positive and we have neglected the the integrals, larger for the latter, hence the two expressions 

absolute value of the Jacobian ^'^^ equivalent if the integrals are dominated by the same 

In addition to Q we consider the definition without '''PP°''^- ^^lus, to extract from p{f) the correct result for 

the theta-functions, which we denote by ptot(/), which ^(/) ^« ^^o^^ld ^e able to isolate the contributions from 

counts the total number of TAP solutions. The effect of the common support. For a generic value of N this could 

theta-functions is to eliminate all solutions with at least ^e quite a hard problem. However, in the limit of large N 

one negative eigenvalue, therefore the meaning of ptot (/) ^"^^'^ ^^'^ integrals are evaluated by saddle point methods, 

is not exactly the same as that of p{f) since all solutions ^ ^™Pl^ ^PP^^^^' 

are now counted. There is just one case in which the two In this case p{f) can be evaluated simply considering 

formulations, at least in the limit of our interest, ^ 1, only the stationary points for which all eigenvalues of the 

are indeed equivalent: if for large N the two integral - Hessian are positive and the Plefka's criterion is satisfied, 

with and without theta-functions - are dominated by the disregarding all others. We stress that such constraint is 
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not contained into p{f), so that the functional alone can- marks the dynamical transition ^Hl- The static transi- 

not give the desired result. tion can be seen as the point where the lowest accessible 

In the next Sections we shall illustrate this procedure (metastable) states dominate. The dynamical transition, 

(re)computing the complexity for the spherical pSP-SG on the contrary, takes place at the point where the be- 

model without external field 9^ using both expressions haviour is ruled by higher, highly degenerate, metastable 

G)) and ®. states. 

The TAP functional has been derived in Refs. [TTII^ : 



3 TAP Equations for the spherical pSP-SG ^i^TAp(m) = ^ V j,^ , m,,---m, 

ii,...,ip 

model ^, /i N 

-y -g) 

The spherical pSP-SG model consists of N continuous — [l + {p ~ l)q^ — pq^^^] (9) 

spins (T, interacting via n-body interactions |H|: i at p i ^ i • i • i- -^i 

^ o J i-i where Nq = }_^ni^, and takmg the derivatives with re- 

N 

H{a) = - ^ erf — ^ Ji^ i ■ ■ - ai (8) spect to one obtains the TAP equations. 

i=l l<ii<-<ip<N 

The structure of the solutions is better understood per- 

The couplings are quenched independent Gaussian vari- 

forming the change of variable mi = q^^^fhi fhf = N) 

ables with zero mean and average (( Jii,....ip)^} — pi/ {2NP^^). 

which leads to TAP functional density: 

The scaling with N ensures a well defined thermodynamic 

/2 

limit. Here and in the following ((•••)) denotes disorder av- /tap(<?, -E) — q E 2^'^^"'^ 

erage. The parameter r is a Lagrange multiplier to impose [l + (p ^ l)?^ ~ ?"Z^ ^] (10) 

the global constraint af ^ N on the spins ampli- ^^^^^ ^ ^ -{l/Np\) ^ J^,,...,^pm,, ■■■m,^ is the T ^ 

energy density. In general, i? is a random variable which 

The study of both the static and dynamical properties depends on both the realization of couplings and on the 
shows that in the thermodynamic limit the model presents orientation of the vector m. However all cases with the 
a (static) transition at a temperature T„ between a high ^^^^ ^^lue of E will also have the same free energy, thus 
temperature replica symmetric phase and a low temper- ^^n consider E as given and study the solutions as a 
ature phase with one step of rephca symmetry breaking function of E. The TAP equations then reduce to 9, /tap (<?, 

Despite its simplicity, the spherical pSP-SG model for q .^^jj^i,^}]^ f^an be written- 

p > 2 has an exponentially large number of locally stable , 

(1 — q) q^' ~ = zT (11) 

states which dominate the dynamical behaviour above Tg. 

where 

As a consequence, two-time correlation functions acquire 

1 



a time persistent part at a temperature Td > Tg which 



p-1 



-E±^/W~E^ , (12) 
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Ec = — ■\/2 (p — l)/p. (13) can be easily seen that the condition (|16|l . for which xp < 

0, leads to an unphysical q decreasing with temperature. 



It is easy to understand that for any positive z and tem- 
perature T below 



T,= il~q,)qPj'-'z-' 



(14) 



where q^ ~ {p — 2)/p, there are two solutions of the TAP 
equation (|ll|l . one larger and one smaller than (/a- Anyway, 
possible candidates for physical solutions are only those 
which are local minima of /tap(<Z, E). 

By using the TAP equation Hll|l the second derivative 
of /tap(9, E) with respect to q evaluated on the solutions 
can be expressed as 

p-2] 



d'jTAP {q,E) 



P 

Apq 

X 



1 



(1-9) 



P 



2 ^Ji{p^l)q^'-^ 



p 

APq 



P-2 
P 



(15) 



where Zc = y^2/p{p - 1) and = l3^p/2. The requirement 
of positiveness of the second derivative thus selects the 
solutions 



P-2 , 
q < tor z > Zc 



q > 



P 
p_-2 

P 



for z < Zc 



(16) 
(17) 



By comparing the two expressions in (|15|l we see that the 
condition z < is equivalent to 



= l-fiip- l)qP-' (1 -q)'>0 



(18) 



which also follows from the stability requirement of the 
replica saddle point |S] and of the dynamics ^Hl- This is 
the Plefka's criterion Q for the physical relevance of TAP 
solutions 13* for the spherical pSP-SG model. Indeed it 



4 Complexity of the spherical pSP-SG model: 
standard calculation 

In this Section we report the main steps of the calculation 
of p{f) for the spherical pSP-SG model. Details can be 
foimd in the literature, see e.g. Refs. [211^1 and for the 
specific case of the p-spin spin glass model Refs. [TSII^II^ . 
The starting point is [cfr. lO] 



rl ^ r p + oo 



drui 5{Gi) 



deiA 



HNq-Y.-n^l] 5 [^TAp(m) - Nf] (19) 



with 



(20) 



and 



N 



where 



a'{q)mi rrij 



1 



a{q) = h/i(p- 1)(1 - q)q 



p-2 



(21) 



(22) 



1-g 

a'{q) — da{q)/dq and we have used the short-hand nota- 
tion: 



^1 1 ■ ■ ■ 7 — 1 

(23) 



+jm ( a] m — iqrn^\. 



(27) 
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and similarly in H21|l . The last term of A is of order 0{1/N), lT ^ T ^^^^^ ^ ^^^^ 
and can be neglected for iV ^ oo [see also below] . 

The structure of the minima is given by the couphngs, j _ /■+°° dm dm f Mg^ ^ (ifh)'^ 
therefore p{f) (and so p{f)) may change from sample to 
sample. Thus, in principle, to have a well defined com- 
plexity we should introduce replicas to compute (lnp(/)) a-nd A — 2/i(p — l)qP ^. 

P] . However it can be shown [HEI that for this model, in The function G.^^ (g) comes from the average of the de- 
absence of a magnetic field, the annealed average ln(p(/)) terminant of A which can be computed either using Grass- 
is exact, so we can just average (O over the disorder. mann variables l3| or introducing replicas pifTSj. The form 

To perform the average over the couplings it is con- depends on the sign of a;p [eq. ^ gO]: 
venient to use the integral representation of the delta- 

G,p(<z) = -ln(l-(?), fora;p>0 (28) 

function to exponentiate its argument. This introduces 

additional parameters which are usually denoted by /, ^^'^ is B ^ solution always adopted in standard calcu- 

q and m m conjugated to /, q and G, and the ad- l^^io^" ^^.g. in and 

ditional variable A coming from Hubbard-Stratonovich n ~ ^ 

" \{l-qf _ 

transformation. The calculation can be further simplified \ 

+ ln-(l-q), forxp<0. (29) 

by substituting ^ Jii,....ip ™ii ' ' ■ ™jp from equation H2U|I 

in fTAp(m) [eq. ®] and noticing that the error involved The two expressions coincide for xp = 0, i.e. for A/2 = 
in disorder-averaging the determinant of A separately ac- 1/(1-9)^- Details of the calculation can be found in the 
counts for changing A of terms of order 0(l/iV) and hence Appendix A. 

negligible as ^ cxo p|[TH| . Integration over m, m, q, A can be done by the saddle 



Performing the averages over the couplings results in Point method 21 , which turns out to be exact for the 

integrals over m,m, A being Gaussian, while the integral 
over / can be easily performed giving a delta-function. 
The final results is then 



m)) = c i df dq dq dAe^'' (24) 



where c is a constant and 



S = if3f[f ^ fiq)] +iqq-A{l-q) 

-jA'' +\nI + G,^{q) (25) 

with f{q) the TAP density functional /tap(»tT') evaluated 
on the solution of the TAP equation H2Q|I : 



(/5(/)) ^c' / d<z4/(g)-/]e^^(«) 
Jo 

^e^^('?*)^ (30) 



where 



/(g) ^-^{i-qn-j{v-m- <i)<f-' m = g^m + ^ + ^ in? - ^ Hm"^') 



A. Crisanti et al.: The Complexity of the Spherical p-spin spin glass model, revisited 



p-1 
2iipqP-'^ 
1 



-fi{l-q)qP-^ 



(31) 



-0.46 



-0.47 



and q* — q*{J) is the solution of 

f{D = f 



(32) 



which gives the largest value of S{q) |22j . 

The simplest way of studying the solution is using q 
as a free parameter to scan all values of /. This is what 
is done, for example, in Ref. [5] where the result (|31|I - H32|I 
was first derived. 

The solutions of eq. H32|) can be found using the re- 
sults of Section 121 The free energy f{q) as function of q 
for all stationary points of /tap(9, E) is shown in Figure^ 
for p — A and temperature T between the static transition 
temperature Tg and the dynamical transition temperature 
Td- Other values of p or T in this range lead to a qualita- 
tively same picture. The corresponding S{q*) as function 
of / is shown in Figure [21 

Strictly speaking to evaluate {p{f )) we should take for 
each value of / the largest value of S , and compute the 
sign of the neglected coefficient in (|30|l . However, one is ac- 
tually interested into the number of metastable states, so 
in all calculations done so far all solutions with Xp < are 
cut out, the "famous" B — solution. Even if not explic- 
itly stated, this is in the spirit of the procedure described 
in Section [3 We stress, however, that if the procedure is 
the same the motivations are not. Indeed the B solu- 
tion can also describe minima of the TAP functional but 



f(q) 



-0.48 



-0.49 



1 

- Xp>0 


1 1 1 1 1 
Xp<0 


1 

\ Xp>0 _ 


— /' 
/ 




\ — 
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Xp<0,q>q^ 

Xp<0,q<q__ 
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0.2 


0.4 0.6 


0.8 



Fig. 1. Free energy density f{q) as function of q for p = 4 and 
temperature T = 0.51 between the static transition tempera- 
ture Ts = 0.5030... and the dynamical transition temperature 
To = 0.5443.... Thicker lines correspond to solutions for which 
the Plefka's criterion is satisfied, while full lines correspond 
to solution for which 1151 is positive, i.e., to local minima of 
/TAp(g,-E'). 

such configurations violate the Plefka criterion (xp < 0), 
thus leading to a non-physical linear susceptibility. 

If the Xp < solutions are disregarded, we are left with 
the curves shown in FigureOcorresponding to the solution 
of the TAP equations with q < (p — 2)/p (dashed line) and 
q > {p ^2)/p (full line). Again, if no other information is 
added, for each / the largest value must be selected to 
evaluate S{f)- This means that there is a region of free 
energies where S{f) is dominated by solutions with q < 
(p — 2) /p. But these are not local minima of /tap, see (|15|) . 
and hence for these free energies S{f) does not give the 
desired result. This clearly shows that the condition Xp > 
alone does not guarantees that only physical states are 
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0.05 



0.04 



E(q*) 



0.03 



0.02 



0.01 



_ Xp>0,q>q^ 

■- ''p>o.q<qa 

--- Xp<0, q>q^ 

— X <0, q<q 
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f 



Fig. 2. S{q*) as as function of / for p = 4 and tem- 
perature T = 0.51 between the static transition tempera- 
ture Ts — 0.5030... and the dynamical transition temperature 
Td = 0.5443.... Thicker lines correspond to solutions for which 
the Plefka's criterion is satisfied, while full lines correspond 
to solution for which H15|l is positive, i.e., to local minima of 
/TAp(g, E). Only values of / for which E > are reported. 
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Fig. 3. Same as Figure |5| where only solutions satisfying the 
Plefka's criterion are displayed. The full line is the result found 
in Ref. 0. 



counted. To find the correct answer additional information 
on the solutions must be added. 

For the spherical pSP-SG this information is easily 
obtained. Indeed, the analysis of the TAP solutions in 
Section 13 shows that only solutions with a;p > and 
q > [p ~2)/p do correspond to metastable states, so only 
the full line in Figure 13 must be considered. This leads to 
the result first derived in Ref. 0. 

5 Complexity of the spherical pSP-SG model: 
Hessian Eigenvalues 

In the previous Section we have revised step-by-step the 
standard calculation of the complexity for the spherical 
pSP-SG model showing which additional information, not 
included into the definition of p(f), must be added to yield 
the correct answer. In this Section we show that that ad- 
ditional information is exactly the thcta-functions needed 
to transform p{f) into p{f)^ see eqs. © and 

To prove the equivalence we must compute the eigen- 
values \i of the Hessian matrix jSJ that for the spherical 
pSP-SG model is given by if^ . 

The eigenvalues are solutions of the equations 

^ A,j = a{q) C, + 2q a' {q) m» ^ ^ rhj 



iSq^P-^^/^ 



A 6 



(33) 



There are two classes of eigenvectors longitudinal 
and transversal. 
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5.1 Longitudinal Eigenvector 

The longitudinal eigenvector is given by: 

and hence satisfies the equation: 



(34) 



[a{q) + 2qa'{q)\ rhi - 



(p-2)! 
At, rhi 



(35) 



Since the Hessian must be evaluated on the solution of the 
TAP equation we can use (|2()|l to write 

^^(p-2)/2 J2 rhP-^ = a{q) {p ~ 1)! m, (36) 
j 

which inserted into H35(l leads to 



XL = 2qa'iq)-ip~2)a{q) 



p-2 



- A*(P - 1)9' 



.p-2 



(37) 



The longitudinal eigenvalue is therefore, apart from posi- 
tive multiplicative coefficients, equal to 9^/tap('Z, E) eval- 
uated in Section O [eq. (|15l) ]. The different coefficients 
come from the derivative being taken with respect to q 
or to mi. The longitudinal eigenvalue has degeneracy 1. 

Note that the term of 0(1/A^) in (|21(l yields a contri- 
bution of 0(1) for longitudinal eigenvectors and cannot 
be neglected [T^. 

5.2 Transversal Eigenvectors 

Transversal eigenvectors satisfy the orthogonality condi- 
tions: 

5^^zTOz = (38) 



and hence span a space of dimension — 1. The eigenval- 
ues equation for transversal eigenvectors can be written 



as 



(39) 



fel<--</Cp-2 



(40) 



For large values of N Jij is a symmetric random matrix 
whose elements are independent Gaussian variables with 
zero average and variance: 



^Iip - l)qP~^ 



N 



(41) 



Therefore for N ^ oo the spectrum of Jij is given by the 
Wigner's semicircular law: |23ll29j . 

p(At) 



27r^(p - l)qP-^ 
X^4ii{p - l)qP-^ - [\t ~ a{q)]^ (42) 

This gives a spectrum at the leading order in N, display- 
ing a non negative support. Since it can be shown that the 
tails of this distribution go to zero exponentially with N 
PS| we can safely exclude negative eigenvalues. The ther- 
modynamic limit transversal fluctuations are, thus, always 
stable, regardless of the sign of xp and the whole stabil- 
ity depends on the longitudinal eigenvalue. Note, however, 
that the A^ — 1 transversal eigenvalues dominate the cal- 
culation of 

det A = exp (Tr In A) for A^ oo (43) 

and any information from the longitudinal eigenvalue is 
washed out when computing p{f). 
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In conclusion we see that the procedure described in 
Section |2 of selecting the saddle point solutions of p{f) 
according to their physical relevance obtained from the 
(independent) analysis of the TAP equations produces the 
correct result for the complexity. 

6 Conclusion and Discussion 

The study of the complex behaviour of glassy systems in 
terms of the topological properties of the energy or free- 
energy surfaces has recently put new interest into the cal- 
culation of the number of metastable states, also called 
complexity or configurational entropy, in mean-field spin 
glass models. In this context classical calculations done 
for the SK and other disordered spin models have been 
reconsidered, extended and also criticised OEIIZIEZI- 

Motivated by these controversies in this paper we have 
reviewed and examined further some questions concern- 
ing the calculation of the complexity of disordered sys- 
tems. Particular care has been taken to distinguish be- 
tween what we would like to compute, and what we 
are able to compute, p{f). We have also discussed how 
information on p{f) can be extracted, at least in the ther- 
modynamic limit. 

The general approach has been illustrated using the 
spherical pSP-SG model, showing the correctness of the re- 
duction procedure. As by-product we have explicitly shown 
that the Plefka criterion separates all solutions of the TAP 
equations into two classes, both containing local minima 
as well as saddles. However, only local minima which sat- 
isfy the Plefka criterion do represent physical states. The 



Spherical p-spin spin glass model, revisited 

Plefka's criterion is indeed a necessary, but not sufficient, 
condition for physical states, and hence cannot be used 
alone for the reduction procedure but the requirement of 
local stability must be added. 

In Ref. 5^ the classical calculation of Bray and Moore 
[2] and DeDominicis and Young 3 for the SK model have 
been critically revised. The main criticism steams from 
the observation that p{f) (identified with p{f) in those 
papers) can be written as saddle point calculation over a 
functional which posses a supersymmetry between com- 
muting and anti-commuting variables used to express the 
Jacobian in {Tj) 26 . The classical solution breaks this sym- 
metry and in Ref. a different, supersymmetric solution 
was proposed. 

In a separate paper J2| we have performed a carefully 
analysis of both the classical and the new supersymmetric 
solutions. The outcome is that both solutions have some 
limitations. For example, the neglected prefactor could be 
exponentially small in N for the classical solution chang- 
ing its prediction |15ll7| . On the other hand the supersym- 
metric solution has a negative Xp, while it is positive for 
the classical solution. This point was recently used [H] to 
justify the correctness of the classical solution. However, 
as we have explicitly shown here for the spherical pSP-SG 
model, this condition is not a sufficient condition. In order 
to prove the correctness of the solution one should prove 
that it corresponds to a physical (stable) state. A rather 
difficult problem already at the annealed (replica sym- 
metric) level used in these calculations. Moreover, since 
it is known that a (marginally) stable solution for the 
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SK model requires an infinite-number of replica symmetry A Appendix 
breakings, the requirement of local stability may partially 
or totally wash-out the results from the annealed approx- 
imation. A complete calculation of the complexity of the 
SK model must include full replica-symmetry breaking. 
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Here we calculate (detA) ioi N ^ oo using the identity 
[see, e.g., Ref. ^] 



detA = Jf[ dri, dr]+ exp | ^ 



(46) 



making not only the calculation but also the analysis of , , + , . , • n • 

where r/i and rj^ are anti-commutmg (Grassmann) vari- 

the saddle points more difficult |19ll27ll3Uj . Supersymme- 
try requirements introduce partial simplifications, since 



ables. From eq. (|21|) we have 



they lead to a connection between the complexity and the 
replica calculation. However, we stress that supersymme- 
try it is not an a-priori requirement for the complexity. 
Indeed, while p{f) is supersymmetric, the reduction pro- 
cedure needed to go from p{f) to p{f) may destroy the 
supersymmetry, so that supersymmetry must be proved 
case by case. It turns out that in the case of the pSP-SG 
model discussed here the supersymmetry is conserved |10J, 
a property which can be associated with the fact that all 
metastable states are uniquely identified by its energy at 
T — and hence its number is conserved |17) . 

Indeed using the TAP equation (lllll S{q) for the spher- 
ical pPS-SG can be rewritten as: 

2 



for xp > and 

+- 



2 — p pz 

In 

p 2 



P - 1 2 
Z 

P 



p'^z^ 



In 



p{p - 1) 



1 



1 



p^p-lfz^ 



(45) 



p(ji — l)z^ 

for xp < 0. Such property is in agreement with the ar- 
gument presented in Sec. I.D of Ref. [7] where the equiv- 
alence between supersymmetry and stability of the TAP 
equations under external field perturbation is shown. 



(detA) = J Y]_d-Vi d-q^ exp ( a{q) ^ 



1=1 

X n 

il<---<ip 



exp 



j n dri, dr]+ exp (^a{q) ^ 

M pKp-^V- 



'7r(ip) 



i=l 



ji<-<« 

X 



2Arp-i (p-2)!2 



(47) 



where J2tt ^ ^^^r all permutations of p different 
integers ii, . . . ,ip. When the products are expanded only 
terms which contain pairs of rjfrji with the same index 
survive. Since since there are (p — 2)! terms with the same 
(44) pairs of Grassmann variables, we end up with 

(detA) = JYI dVidrit exp ^a{q)Y^ri+ri}j 
M p!(p-l)! 



i=i 
X exp 



27VP-1 {p~2)l 



N 



Y]^ drii drj. 



Ji<-< 

+ 



exp 



2N ^ 



i 



(48) 
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